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I.  INTRODUCTION 


The  removal  of  heat  from  electronic  equipment  to  assure 
reliable  operation  has  been  a  problem  plaguing  the  electronic 
industry  for  many  years.  As  a  result,  many  researchers  have 
attempted  to  develop  specific  correlations  to  enable  designers  to 
attempt  to  predict  possible  solutions  to  this  problem. 

This  research  concerns  air  flowing  across  cylindrical 
elements  mounted  in  an  array  on  one  wall  of  a  vertical  adiabatic 
channel.  This  configuration  is  much  the  same  as  an  array  of 
transistors  in  the  T05  configuration  arranged  on  stacks  of 
parallel  circuit  boards  in  various  types  of  electronic  equipment. 

A.  THEORY  OF  VERTICAL  CHANNELS  FORMED  BY  CIRCUIT  BOARDS 
Many  electronic  systems  consist  partially  of  arrays  of 
printed  circuit  boards  supporting  a  fairly  dense  array  of 
electronic  elements  on  one  side  as  shown  in  Figure  1.1.  These 
circuit  hoards  can  be  arranged  so  that  the  side  with  the  elements 
( r  n  n  t  sidel  faces  the  relatively  smooth  backside  of  another. 

^ h o n  ^  series  of  circuit  boards  are  arranged  in  such  a  fashion 
Q  ■■i'j'jro  1.2),  long  vertical  channels  in  which  cooling  air  can 
-p  s  , - s 1  < p i  by  forced  convection  can  result.  Durinq  ooeration, 
thp  pip-ip^ts  such  as  flatoacks,  dips  resistors  and  transistors 
generate  neat,  the  majority  of  which  is  dissipated  to  the 
sur round i ngs  with  little  transmitted  back  through  the  circuit 
board.  As  a  result,  it  apoears  as  if  all  of  the  heat  is 
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transmitted  by  one  side  of  the  channel  or  one  side  of  each 
circuit  board  to  the  air  stream  or  air  layer  between  the  boards. 

Complexity  is  added  to  the  system  when  elements  serve 
different  purposes  and,  as  a  result,  possess  different  sizes, 
shapes  and  heat  dissipating  characteristics. 

The  process  of  heat  removal  from  electronic  systems  creates  a 
very  comlex  design  problem  governed  by  such  factors  as  circuit 
board/element  arrangement,  spacing,  heating  characteristics  and 
operating  environment.  Moreover,  the  ultimate  goal  of  any 
thermal  analysis  and  the  prediction  of  operating  temperatures  is 
the  assessment  of  reliability.  Reliability  is  the  probability 
that  an  electronic  component/system  will  operate  properly  without 
failure  for  its  intended  lifetime  under  its  environmental 
operating  conditions.  [Ref.  1] 

For  many  years,  there  have  been  numerous  attempts  to  derive 
techniques  to  aid  designers  in  this  field.  Unfortunately,  the 
methods  developed  are  often  too  complex,  inadequate  or  restricted 
to  specific  conditions. 

Early  work  with  natural  convection  between  parallel  arrays  of 
vertical  circuit  boards  was  performed  by  Elenbaas  [Ref.  2], 

Bodoia  and  Osterle  [Ref.  3]  and  Engel  and  Muller  [Ref.  4].  In 
1972,  Aung  et  al .  [Ref.  5]  studied  the  effects  of  circuit  board 
supports  and  elements  in  the  field  of  flow  for  wide  channel 
spacinqs.  Aung  went  on  to  investigate  uniform  heat  flux  on  two 
channel  walls  in  which  he  subsequently  studied  the  effects  of 
baffles  and  circuit  board  staggering.  [Ref.  6]  There  is 
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considerable  literature  describing  the  efforts  that  have  been 
made  to  attempt  to  model  circuit  channel  cooling  by  natural 
convection  and  forced  convection. 

For  example,  Kelleher  utilized  a  thermal  network  approach  to 
derive  a  more  repr esentative  model  of  forced  convection  for  o.ne 
specific  qeometry  with  equal  heat  releases  [Ref.  7],  Other 
experimental  or  semi -expe rimental  approaches  were  performed  by 
Tamura  [Ref.  8],  Finch  [Ref.  9],  Laermer  [Ref.  10],  Marto  [Ref. 
11]  and  Parkes  and  Preston  [Ref.  12].  Most  of  this  work  is 
restricted  to  specific  conditions. 

B.  THEORY  OF  SURFACE  ROUGHNESS  GEOMETRIES 

Surface  roughnesses  have  been  classified  by  several  different 
techniques  of  which  the  most  commonly  used  are  the  k  and  d  type 
roughness  classifications  consisting  of  2  and  3  dimension  classes 
which  are  further  divided  into  artificial  and  commercial 
roughness  subclasses.  Although  many  researchers  allude  to  these 
classifications,  only  a  limited  number  actually  use  this 
terminology  in  agreement.  It  is  further  noted  that  in  the  area 
of  heat-transfer,  the  k-  and  d-  type  roughnesses  seem  irrelevant 
because  the  same  heat  transfer  roughness  parameter  (g+) 
correlation  is  proposed  for  all  types  of  roughnesses. 

C.  THEORY  OF  HEAT  TRANSFER  FROM  ROUGH  SURFACES 

Earlier  researchers  in  this  field  seem  to  agree  that  the 
heat-transfer  roughness  parameter  (g+),  commonly  referred  to  as 
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the  inverse  of  the  sublayer  Stanton  number,  takes  the  functional 
f  o  rm  : 

P  P 

3+  =  1 / S ts 1  ay e r  =  C(e  )  (Pr) 

where  e+  is  the  roughness  Reynolds  number  based  on  an  arbitrary 
characteristic  height  of  the  roughness  (e).  The  different  values 
of  p  and  q  fall  into  two  categories: 

1.  The  assumption  of  similarity  near  the  wall  which  results 
in  value  of  p  =  .2  and  q  =  .44  or  for  the  eddy-di ffusivi ty 
models  p  =  .25  and  q  =  .66  [Refs.  13  and  14]. 

2.  The  assumption  of  a  breakdown  in  close  to  the  wall 
similarity  results  in  values  of  p  =  .45  and  q  =  .57  or 
p  =  .5  and  q  =  .66  for  the  eddy  diffusivity  models 
[Refs.  1 5  and  14] . 

With  regard  to  this  study,  the  foregoing  is  only  of  academic 
interest  because  the  geometry  here  is  defined  as  an  in-line  array 
of  cylindrical  elements. 


1 1 .  EQUIPMENT  DESIGN  AND  INSTRUMENTATION 


The  equipment  utilized  in  this  research  was  either 
manufactured  by  the  Engineering  Machine  Shop  or  supplied  by  the 
Naval  Engineering  Department  of  the  Naval  Postgraduate  School., 
Monterey,  California. 

A.  WIND  TUNNEL 

The  wind  tunnel  (Figure  2.1)  is  constructed  of  1/8" 
plexiglass  with  the  exception  of  two  1/4"  access  plates  and 
various  reinforcement  members  mounted  on  the  exterior  so  as  to 
ensure  smooth  air  flow  on  the  tunnel  interior. 

The  test  section  of  the  tunnel  has  a  12"  by  4"  cross  section 
and  a  height  of  18"  with  access  covers  on  both  faces  to 
facilitate  instrumentation  and  support  for  the  test  surfaces. 

The  tunnel  intake  is  belled  with  an  18"  radius  to  produce  an 
orifice  that  has  a  40"  by  12"  cross  section.  An  adjustable  leg 
is  attached  at  each  corner  to  serve  as  a  support  for  the  tunnel, 
to  level  the  tunnel  and  to  provide  an  8"  clearance  off  the  floor. 

The  upper  end  of  the  tunnel  is  reduced  to  a  section  of  2"  PVC 
iuctinq  which  is  connected  to  a  550  cfm  turbo  fan.  Air  flow 
through  the  tunnel  is  controlled  by  a  series  of  valves  connected 
in  line  with  the  fan. 

A  1/8"  plexiglass  splitter  olate  was  mounted  1/2"  off  of  the 
access  plate  facing  the  test  surface  in  order  to  reduce  the 
effects  of  the  boundary  layer  formed  in  the  tunnel.  This 
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distance  was  pre-de termi ned  through  the  use  of  computer  software 
[Ref.  16]  and  the  results  are  displayed  in  Figure  2.2.  A  wooden 
frame  box  (4'x4'x  8')  was  constructed  around  the  wind  tunnel  to 
elliminate  the  effects  of  incident  radiation  on  the  cylinders 
from  the  laboratory  environment  and  the  box  also  served  to  reduce 
any  transient  air  currents  created  by  movement  in  the  room  (i.e., 
doors  opening  and  closing). 

B.  TEST  SURFACE 

The  design  of  the  test  surface  used  in  this  research  was 
derived  from  the  previous  work  of  Arvizu  and  Moffat  [Ref.  17]  and 
Ortega  [Ref.  18].  It  was  constructed  of  a  12"  x  12.5"  x  1/8" 
sheet  of  plexiglass  (for  support)  with  a  12"  x  12.5"  x  3/8"  sheet 
of  balsa  wood  (thermal  boundary)  cemented  to  it.  The  leading 
edge  of  this  (balsa  wood- pi  exi  g  1  a  s  s )  splitter  plate  was  cut  at  a 
35°  anqle  to  form  a  sharp  leading  edge  to  reduce  the  effects  of 
leading  edge  separation.  The  surface  (balsa  wood)  of  the 
splitter  plate  was  then  sanded  and  painted  three  times  with  white 
enamel  paint  and  sanded  a  fourth  time  to  ensure  a  smooth  test 
surface.  Small  (1/16"  diameter)  holes  were  drilled  through  the 
plate  to  permit  the  installation  of  resistor  leads  and 
thermocouples  into  the  cylindrical  elements  mounted  on  the  balsa 
wood.  Three  mounting  screws  were  attached  to  the  plexiglass  side 
of  the  splitter  plate  and  to  the  access  plate  to  facilitate  the 
positioning  of  the  test  surface.  A  tail  fin  mounted  to  the  tail 
end  of  the  splitter  p1ate  was  used  to  regulate  the  air  flow  on 
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both  side  of  the  splitter  plate.  A  fourth  adjusting  screw  was 
attached  to  this  fin  and  through  the  access  plate 
(Figure  2.3). 

C.  CYLINDRICAL  ELEMENTS 

The  cylindrical  elements  used  in  this  research  were 
manufactured  from  3/4”  round  aluminum  stock  and  cut  to  a  height 
of  1.0  cm.  A  7/16"  diameter  hole  was  cut  to  a  depth  of  1/4"  in 
the  center  of  the  element  and  two  additional  holes  (1/8" 
diameter,  1/4  deep)  were  cut  opposite  each  other  to  permit  the 
installation  of  the  56  ohm,  1/2  watt  resistors.  Another  hole 
(1/32"  diameter,  3/8"  deep)  was  drilled  to  position  the 
thermocouple  (Figure  2.4). 

The  resistors  were  epoxied  into  forty  two  of  the  cylinders 
and  six  more  were  cut  in  half  so  that  an  array  of  7x6  cylinders 
could  be  mounted  evenly  on  the  splitter  plate  and  the  12  half 
cylinders  could  be  arranged  along  the  edges  of  the  plates  (Figure 
2.5)  . 

The  resistor  (56  ohm,  1/2  watt),  prior  to  insertion  into  the 
cylinders,  had  the  leads  cut  and  smaller  leads  (2841/1,  *30  solid 
spc.)  attached  so  as  to  reduce  heat  losses  by  conduction  through 
them. 

The  cylinders  were  mounted  on  the  test  surface  in  an  in-line 
array  with  a  separation  of  one  cylinder  diameter  between 
cylinders.  This  was  done  both  vertically  and  horizontally. 


Flow  Adjustment 


Figure  2.3  Sketch  of  Test  Surface  Configuration. 


Top  Row:  Assorted  Views  of  Cylinders 
Bottom  Row:  T05  Transistor  and  a  nickel 


Figure  2.4  Photograph  of  Cylinders 


23 


INSTRUMENTATION 


1 .  Temperature  Measurements 

Type  T,  copper-constantan  thermocouples  were  utilized  for 

the  determination  of  the  temperature  of  the  heated  elements  and 
the  inlet  and  outlet  temperatures  of  the  wind  tunnel  near  the 
test  surface.  The  thermocouples  were  connected  to  seven  scanner 
boards  in  an  autodata  nine  and  calibrated  to  an  accuracy  of 
j+0.2°F  (Appendix  A.)  The  thermocouples  were  then  routed  through 
an  access  plate  in  the  tunnel  through  the  splitter  plate  and  into 
each  of  the  el  erne nts . 

2 .  Velocity  Measurements 

The  velocity  measurements  were  obtained  by  the  use  of  a 
micro  manometer  and  a  pitot-static  tube.  Changes  in  the  dynamic 
and  static  pressure  were  measured  by  the  micro  manometer  and 
using  Bernoulli's  equation,  velocity  was  determined.  The  pitot- 
static  tube  was  placed  between  the  splitter  plate  and  the  test 
surface  splitter  plate  and  used  to  determine  the  approach 
velocity  between  them. 

A  TSI  (Model  1050/1050AA)  hot  wire  anemometer  was  used  to 
determine  the  velocity  profile  prior  to  entry  into  the  test  area. 
The  hot  wire  anemometer  was  calibrated  with  the  micro  manometer 
and  used  only  to  ensure  that  the  approach  velocity  was  uniform. 

Accuracy  of  the  mirco  manometer  was  determined  to  be 
+.0005  inches  of  water  and  resulted  in  an  error  of  +0.2  ft/sec  or 


+7cm/s.  because  the  micro  manometer  was  utilized  for  the 


calibration  of  the  hot  wire  anemometer,  more  accuracy  might 
result  if  the  micro  manometer  was  used  to  determine  the  velocity 
between  the  splitter  plates. 

3 .  Power  Measurements  and  Control 

A  LAMBDA  (60  volt,  10  amp)  Regulated  Power  Supply  (Model 
LK  345AFM)  was  used  to  power  to  the  56  ohm,  1/2  watt  resistors 
mounted  in  the  forty  two  cylindrical  elements. 

A  switch  board  was  set  up  to  allow  power  to  be  delivered 
to  individual  elements  in  either  of  two  columns,  the  third  or  the 
fourth  columns.  It  was  also  set  up  in  such  a  manner  that  each  or 
all  rows  could  be  energized. 

A  Hewlett  Packard  digital  multi  me  ter  (346 6A)  was  used  to 
measure  the  voltage  drop  across  the  resistors  and  a  Fluke  8060A 
multimeter  was  used  to  measure  the  current  through  the  circuit. 
The  power  delivered  is  the  product  of  current  and  voltage. 

4 .  Data  Acquisition 

The  Autodata  Nine  mentioned  above  is  a  self  contained 
recordi nq  device  with  up  to  100  channels  available  for  both 
temperature  and  voltage  measurements.  During  this  research,  71 
channels  were  utilized,  one  for  voltage  measurement  and  the  rest 
for  temperature  measurements  in  degrees  Fahrenheit.  (Figure  2.6) 

E.  EXPERIMENTAL  PROCEDURE 

For  this  research,  data  was  recorded  at  three  different 
channel  heiqht  to  cylinder  height  ratios  (1.0,  2.3  and  4.6)  for 
two  different  approach  velocities  (1.5  m/s  and  3.0  m/s).  The 


data  recorded  consisted  of  the  voltage  supplied  to  the  circuit, 
the  current  through  the  circuit  and  the  temperatures  of  the  test 
elements  as  well  as  the  inlet  and  exhaust  temperatures  of  the 
airstream.  For  each  channel  height  and  approach  velocity,  data 
was  gathered  after  a  steady  state  condition  was  reached,  for 
each  element  of  a  column  energized  one  at  a  time  and  for  the 
first  four  rows  and  for  all  six  rows  energized  all  at  once. 


III.  DATA  REDUCTION 


A.  DATA  TAKEN 

Inlet  Free  Stream  Temperature,  T^  (°F) 

Cylinder  Surface  temperature,  Ts,0  f°F) 
Atmospheric  Pressure,  Patm  (inches  Hg.) 
Pitot-Static  Pressure  Difference,  AP  (inches  H 2 0 ) 
Power  Supply  Voltage,  V  (volts) 

Current  Drawn,  I  (amps) 

T s , i  ( °F  ) 


B.  SOUGHT  AFTER  PARAMETERS 

1 .  Dimensionless  Temperature  Ratio 

The  Dimensionless  Tempe ra ture Rati  0 ,  0,  is  obtained  as  a 
function  of  cylinder  location,  channel  configuration  and 
freestream  velocity.  Here  9  is  defined  as 


0  a  Ts,i  '  T- 

i  T  ,  -  T 
s ,  1  °° 

where  Ts,i  is  the  surface  temperature  obtained  by  the  first 
passive  cylinder  directly  after  the  cylinder  heated  by  a 
resistor.  Ts,i  is  the  surface  temperature  obtained  by  individual 
passive  cylinders  after  and  in  the  same  column  as  the  cylinder 
heated  by  the  resistor  where  i  is  the  row  after  the  heated 
cy 1 i nde  r .  Note  that 


a 

1 


1.0  by  definition 


2.  Heat- transfer  Coefficient 


The  Heat- trans fer  Coefficient,  h(W/m2°C),  is  found  and 
is  a  function  of  the  cylinder  location,  channel  configuration  and 
freestream  air  velocity. 

The  heat  transfer  coefficient  was  determined  as  follows: 

q  ( W)  -  Q  ( W ) 

h=  input  losses 

- A  TTs'.o  -  T»  ) - 

Ts,0  and  T^  were  recorded  by  the  autodata  nine  and 
corrected  by  the  calibration  data  in  Aopendix  A. 

A,  is  the  surface  area  of  the  cylinder  is  calculated  for 
channel/cylinder  height  ratios  of  1.0,  2.3  and  4.6.  If  the 
channel /cyl i nder  height  ratio  is  equal  to  1.0,  A  =  t d L  and  if  the 
c hannel /cy 1 i nde r  height  ratio  is  greater  than  1.0,  A  =  -[dL  + 
d2/4]  . 

The  value  of  Qiosses  in  watts  is  derived  from  a  parabolic 
temperature  dependent  approximation  utilized  by  Ortega  [Ref.  18] 
in  his  researach  with  aluminum  cubes.  This  approximation  is 
assumed  to  be  the  same  for  cylinders  and  is  used  here. 

^losses  (w)  =  "0,62927  x  10’3  +  0.5059  x  10"^(Ts,o-Tx) 

-  0.023  695  x  1  0-5  (Ts,0-Tj2 

Ortega  also  determined  the  radiation  losses  to  be  an  order  of 
magnitude  smaller  than  those  given  above,  and  radiation  losses 
were  therefore  disregarded. 


QinputfW)  is  the  product  of  power  supply  voltage  and 


current  drawn. 


Q input  =  V( volts)  x  I(amps) 

3 .  Reynolds  Number 

The  Reynolds  number  is  developed  as  a  function  of  channel 
configuration  and  free  stream  air  velocity.  In  this  research, 
the  Reynolds  number  was  calculated  using  the  freestream  air 
velocity,  channel  height  and  viscosity  as  follows: 


Re  =  channel  height  x  velocity 
H  vi sc  os  i  ty 


The  value  of  the  viscosity  is  obtained  by  a  linear 
interpolation  of  data  from  tables  of  viscosity  and  air 
temperature. 

4  .  Nusselt  Number 

The  Nusselt  number  is  prescribed  as  a  function  of  channel 
configuration.  Here,  the  Nusselt  number  was  calculated  using 
channel  heiqht,  heat- trans f er  coefficient  and  the  thermal 
conductivity  of  the  free  stream  air  as  follows: 


N u  =  channel  height  x  heat-transfer  coefficient 
H  thermal  c  o  n  d  u  c  t i v i ty 


The  value  of  the  thermal  conductivity  was  obtained  by  a 
linear  i  n te rpo  1  a t i on  of  data  from  tables  of  thermal  conductivity 


and  air  temperature. 


c.  typical  set  of  data 

Table  I  presents  a  typical  set  of  data  obtained.  These  data 
will  be  used  here  to  illustrate  the  data  reduction  procedure. 

TABLE  I 
DATA  SET 

patm  =  29.736  inches  Hg. 

LP  =  0.0051  inches  H2O 
Channel/Cylinder  Height  Ratio  =  4.6 
Cylinder  Heated  -  row  1,  column  3 
Cylinder  Surface  Temperature,  Ts,q  =  206. 1°F 


Vol tage 

Su  ppl i ed 

* 

14.380 

volts. 

Current 

D  r  awn 

= 

0.225 

amps. 

Ts»l 

s 

80 . 3°  F 

T  s  >  2 

= 

7  4 . 0°  F 

T  s » 3 

= 

7  2 . 2  0  F 

Ts.4 

= 

71 .5°F 

T  s  >  5 

= 

70 . 8° F 

= 

68. 6°F 

D.  CALCULATIONS 

1 .  Dimensionless  Temperature  Ratio, 

The  temperatures  for  the  cylinder  surface  were  recorded 
by  the  autod.ta  nine  and  corrected  by  use  of  the  thermocouole 
calibration  data  presented  in  Appendix  A. 

T  s ,  -j  (corrected)  =  m  x  Ts,-j  (recorded)  -  b 
m  Recorded  j)  T  s ,  i  (Corrected) 


Tc.l 


1.0026  (30.3)  -  0.9038 


79.6  F 


Ts,0  =  1.0025  (206. 1°F)  -  0.4745  =  206. 1°F  (corrected) 

T,  =  1.0018  ( 68 . 6 ° F )  -  0.5980  =  68.1°F  (corrected) 

Ts>0  -  T,  =  206.1  -  68.1  =  138. 0°F  =  76.7°C 


s 


Channel/cylinder  height  ratio  =  4.6 
therefore  A  =  -[(0.019m)  (0.01m)  +  (0.O19m)?/4] 
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A  =  0.00088m2 

Olosses  =  0.692  7  x  10“3  +  0.5059  x  10-2(TS>0-TCO) 

-  0.23695  x  10-  6(Ts,o-tJ2 

where  Ts,0  -Too  is  in  °C  and  Olosses  1S  in  watts. 

Ts  ,  0-  T  oo  =  76.7  C 
Olosses  =  0.373  watts 

Q input  =  =  14.380  volts  x  0.225  amps  =  3.27  watts 

h  =  3.27  -  0.373  =  42.92  W/m2  °C 

OTD'0'088  (7  6.7) 

3 .  Reynolds  Number 

ReH  =  HU^/  v 

H  =  4. 6  ( .01m)  =  0.046m 
U  =  1  .  5m/  s 

=  38.02  x  lO'6  -  109.8  x  10'9  (27.0  -  T  x  ) 
where  Ty,  is  in  °C 
T  =  2  0 . 0  °  C 

GO 

=  38.02  x  10-6m2/s 

Re^  =  .046m  x  1.5m/s  =  1814.7 

38.02  x  10“6m2/s 

4  .  Nu  ssel t  Number 


Nuh  =  Hh/k 

H  =  ,046m  h  =  42.92  W/m2°C 

k  =  0.0407  -  76.8  x  10“6  (27.0-T  ) 


IV.  RESULTS 


The  experimental  measurements  that  were  obtained  lead  to  two 
different  subject  categories:  thermal  wake  effects  and 
heat-transfer  coefficient.  Table  II  presents  the  range  of  the 
parameters  covered  and  a  listing  of  the  data  is  tabulated  in 
Appendix  C. 


TABLE  II 

PARAMETERS  TESTED 


Thermal  Wake 


Channel  Height/Cylinder  Height  Ratio  1.0  4.6 

Downstream  Row  Number  0,1, 2, 5  0,1, 2, 5 

Velocity  (m/s) _ 1.5,  3.0 _ 1.5,  3.0 


Heat-Transfer  Coefficient 


Channel  Height/Cylinder  Height 

Ratio  1.0 

2.3 

4.6 

Heated  Cylinder  Row  Number 

1,2, 3, 6 

1  ,2,3,6 

1,2, 3, 6 

Number  of  Rows  Heated 

4,6 

4,6 

4,6 

Vel oci ty  (m/s) 

1.5,  3.0 

1.5,  3.0 

1.5,  3.0 

A.  THERMAL  WAKE 

The  thermal  wake  created  by  a  single  heated  cylinder  in  an 
array  is  influenced  by  the  channel  height  to  cylinder  height 
ratio,  the  mass  flow  rate  and  the  position  of  the  cylinder  in 
the  array. 

For  a  channel  height/cylinder  height  ratio  of  1.0,  there,  is 
no  free  stream,  and  the  flow  through  the  array  undergoes  periodic 
accelerations  and  decelerations.  It  therefore  seems  unlikely 
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that  a  channel  heiqht/cyl  i  nder  height  rat1-0  of  4,5  would  possess 
the  same  convective  means  for  determination  of  the  downstream 
wake.  However,  a  plot  of  5  as  a  function  of  downstream  row 
number  (Fioure  4.1  and  4.2)  demonstrates  a  1 / N  dependence  for 
both  velocities  at  both  channel  heights. 

Figure  4.3  shows  an  even  closer  resemblance  to  the  1  / N  curve 
for  the  thermal  wake  created  by  a  heated  cylinder  located  in  a 
row  internal  to  the  matrix  of  cylinders. 

The  following  is  provided  for  the  understanding  of  the 
graphical  curves  found  in  the  legends  of  various  graphs  in 
Chapter  I V : 

Hx.xRxVx.x 

Hx.x  Channel  Height/ Cylinder  Height  Ratio  x.x 

Rx  Row  that  the  heated  cylinder  is  in  x 

V x . x  Approach  Velocity  x.x 

For  example,  H.4.6R1V3.0  denotes  data  gathered  at  a  channel 
heiqht/cyl i nder  height  ratio  of  4.6,  where  the  cylinder  heated 
is  in  the  first  row  and  the  approach  velocity  is  3.0  m/s. 

B.  HEAT-TRANSFER  COEFFICIENT 

In  Figure  4.4  the  heat-transfer  coefficient,  h,  for  a  single 
heated  cylinder  in  the  first  row  is  plotted  against  channel 
height/cylinder  height  ratios  for  approach  velocities  of  1.5  m/s 
and  3.0  m/s.  It  can  be  seen  that  as  the  channel  height 
increases,  the  value  of  h  drops  off  fairly  rapidly  to  an 
asymototic  limit  at  a  channel  height /cylinder  height  ratio  of 
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Figure  4.2  Normalized  Dimensionless  Temperature  vs.  Row 

Position  for  a  Channel  Height/Cylinder  Height  Ratio 
of  4.6  and  Heated  Element  in  the  First  Row. 


Figure  4.4  Heat- trans f er  Coefficient  vs.  Channel  Height/ 
Cylinder  Height  Ratio  With  No  Rows  Upstream 


4.6.  Air  flows  over  the  cylinders  instead  of  through  them  (the 
path  of  least  resistance). 

Figure  4.5  show  a  plot  of  h  against  the  channel  height  to 
cylinder  height  ratio  for  an  internal  cylinder  where  h  drops  off 
a  1  i  ttl  e  more  rapidly  . 

Figure  4.6  plots  h  against  array  position  for  a  channel 
height/cylinder  height  ratio  of  1.0.  Here,  h  increases  with 
position  upstream.  This  correlates  to  the  increased  turbulence 
with  position  downstream. 

Figure  4.7  is  a  plot  of  h  against  array  position  for  a 
channel  height  to  cylinder  height  ratio  of  4.6.  For  this 
configuration  h  appears  to  increase  slightly  and  then  drop  off 
slightly  so  that  the  value  of  h  at  the  sixth  row  is  less  then 
that  of  the  first  row. 

Figure  4.8  plots  Nusselt  number  against  Reynolds  number  for 
all  three  channel  height  to  cylinder  height  ratios.  This 
composition  leads  to  a  correlation  of  the  form  Nu  =  C 
where  C  is  different  for  different  array  densities  as  noted  by 
Arvizu  and  Moffat  [Ref.  17]. 
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ONE  ROW  UPSTREAM 


Figure  4.6  Heat- transfer  Coefficient  vs.  Number  of  Upstream 
Rows  for  a  Channel  Height/Cylinder  Height  Ratio 


Rows  for  a  Channel  Height/Cyl 


Figure  4.8  Nusselt  Number  vs.  Reynold  Number  for  all  Three 
Channel  Height  to  Cylinder  Height  Ratios. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  research,  for  cylindrical  elements  mounted  in  a 
vertical  channel,  the  characteri stic  curves  produced  closely 
resemble  those  for  a  similar  configuration  of  cubical  elements 
which  were  investigated  by  of  Arvizu  and  Moffat  [Ref.  17]. 

Although  there  is  room  for  more  data  to  assist  in  making  firm 
conclusions  about  the  flow  characteri  sties  of  air  flow  around 
cylinders,  the 

The  following  has  been  observed: 

1.  The  1/N  characteri stic  curve  produced  by  the  thermal  wake 
of  a  heated  cylinder  on  an  adiabatic  wall  in  a  channel, 
may  allow  for  fairly  accurate  prediction  of  temperature 
increase  in  downstream  elements  in  an  array  for  various  air 
fl ow  rates. 

2.  Channel  spacings  and  airflow  rates  have  very  little  affect 
on  the  shape  of  the  characteri stic  curve  produced  by  the 
thermal  wake  of  a  heated  cylinder  in  an  array. 

3.  The  position  of  a  heated  cylinder  in  an  array,  on  an 
adiabatic  wall  in  a  channel,  has  a  definite  affect  on  its 
heat  transfer  characteristics.  Elements  requiring  more 
cooling  should  be  placed  closer  to  the  freestream  air  inlet 
This  would  reduce  effects  created  by  the  thermal  wake  of 
other  cylinders. 

4.  Heat  transfer  in  an  array  of  cylindrical  elements  is 
mainly  influenced  by  the  air  flow  between  the  cylinders 
and  slightly  less  by  the  turbulence  in  the  air  flow. 

5.  As  the  channel  height  to  evlinder  height  increases,  the 
amount  of  heat  transfer  that  occurs  decreases  to  an 
asymototic  limit  near  a  channel / c  y 1 i nder  height  ratio  of 
4.6. 

For  continued  work  in  this  field,  it  is  recommended  that,  a 
longer  test  surface  with  more  rows  of  elements  be  utilized  along 
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APPENDIX  A 


THERMOCOUPLE  CALIBRATION 

The  thermocouples  were  all  connected  to  a  wooden  framed 
switchboard  which  was  connected  by  10'  wire  leads  to  seven 
scanner  boards  mounted  in  the  autodata  nine  data  acquisition 
device.  This  combination  was  then  calibrated  as  a  complete  unit 
in  the  calibration  laboratory  against  five  known  temperatures : 

6 6 . 0 2° F ,  1 06. 14°F ,  127.615°F,  151.20°F  and  200.42°F 
Thermocouple  readings  follow: 


1 

44.5 

104.4 

127.9 

151.5 

200.5 

2 

44.5 

104.5 

128.0 

151.5 

200.4 

5 

44.2 

104. S 

128.0 

151.2 

200.4 

4 

44.5 

104. 5 

128.0 

151.5 

200.4 

5 

44.1 

104.4 

127.9 

151.2 

200.4 

4 

44.2 

104.4 

127.9 

151.2 

200.4 

7 

44.1 

104.2 

127.8 

151.1 

200.5 

8 

44.2 

104.4 

127.8 

151.5 

200.4 

9 

44.1 

104.2 

127.7 

151.0 

200.2 

10 

44.2 

104.4 

127.8 

151  .2 

200.4 

n 

44.8 

104.9 

128.5 

151.8 

201.0 

12 

44.8 

107.0 

128.5 

151.8 

200.9 

15 

44.7 

104.8 

128.5 

151.4 

200.8 

14 

44.7 

107.0 

128.4 

151 .8 

201.0 

IS 

44.4 

104.8 

128.5 

151.4 

200.7 

14 

44.7 

104.8 

128.5 

151.7 

200.8 

17 

44. 5 

104.7 

128.0 

151.4 

200.7 

18 

44.7 

104.8 

128.2 

151.4 

200.8 

19 

44.  S 

104.4 

128.0 

151.4 

200.7 

20 

44.2 

104.8 

122.4 

152.4 

199.5 

21 

44.2 

104.4 

128.0 

151.5 

200.4 

22 

44.5 

104. S 

128.0 

151.5 

200.5 

25 

44.1 

104.5 

127.9 

151 .2 

200.4 

24 

44.2 

104.5 

128.0 

151.5 

200.5 

25 

44.0 

105.5 

127.8 

151.2 

200.5 

24 

44.2 

104.5 

127.8 

151.2 

200.4 

27 

44.0 

104.2 

127.1 

151.0 

200.2 

28 

44. 1 

104.2 

127.7 

151.1 

200.2 

29 

44.9 

104.2 

127.7 

151 . 1 

200.5 

50 

44.5 

104.5 

127.8 

151.2 

200.  1 

51 

44.5 

104.4 

127.8 

151.2 

200.4 

52 

44.5 

104.5 

128.0 

151.4 

200 . 4 

55 

44.2 

104.5 

127.9 

151.4 

200.6 

54 

44.5 

104.5 

128.0 

151.4 

200.4 

55 

44.5 

104.5 

127.9 

151.2 

200.4 

54 

44.5 

104.5 

128.0 

151.4 

200.5 

57 

64.5 

104.5 

127.9 

151.5 

200.5 

58 

44.4 

104.5 

128.0 

151.5 

200.4 

59 

44.5 

104.4 

127.9 

151.2 

200.4 

40 

44.4 

106.7 

128.2 

151.4 

200.7 

41 

44.5 

104.4 

128.1 

151.5 

200.7 

42 

64.6 

104.8 

128.4 

151.4 

200.9 

45 

64.4 

106.7 

128.2 

151.5 

200.8 

44 

64.7 

104.8 

128.4 

151.4 

200.8 

45 

64.5 

105.8 

128.2 

151.5 

200.7 

44 

64.6 

106.8 

128.5 

151.6 

200.8 

47 

64.5 

106.7 

128.2 

151.6 

200.7 

48 

64.7 

104.7 

128.2 

151.6 

200.7 

49 

66.4 

104.4 

128.0 

151.5 

200.7 

50 

44.4 

104.5 

128.1 

151.4 

200.5 

51 

44.5 

104.4 

127.9 

151.5 

200.6 

52 

44.4 

104.6 

128.2 

151.4 

200.7 

55 

64.5 

104.6 

128.1 

151.4 

200.4 

54 

64.4 

104.6 

128.2 

151.4 

200.7 

55 

66.4 

106.4 

128.0 

151.4 

200.6 

54 

44.4 

106.6 

128.5 

151.5 

200.7 

57 

64.4 

104.5 

128.0 

151.4 

200.6 

58 

64.4 

106.6 

128.2 

151.4 

200.6 

59 

66.5 

104.6 

128.1 

151.4 

200.6 

40 

46.2 

106.8 

128.0 

151.5 

200.4 

41 

66.2 

104.7 

128.0 

151,4 

200.5 

62 

44.5 

107.1 

128.2 

151.4 

200.6 

45 

44.5 

106.8 

128.0 

151.5 

200.6 

64 

44.5 

106.9 

128.2 

151.4 

200.6 

45 

44.5 

106.8 

128.1 

151.4 

200.6 

64 

44.5 

104.8 

128.1 

151.4 

200.4 

67 

44.4 

104.9 

128.5 

151.5 

200.7 

68 

66.5 

106.8 

128.1 

151.5 

200.7 
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A  least  squares  fit  program  was  developed  to  produce 
thermocouple  corrections  with  the  following  result: 


K 

£ 

r 

i 

l 

l 

i 

o 

r, 

j- 

r 

s 

* 

i 

t 

> 

J 
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"’’actual  *  M  x  ^measured  +  8 


THERMOCOUPLE  l 

Mf* 

1.0017255 

I* 

*0.4258500 

THeFMOCOUtLE  JS 

N* 

1.0028455 

>* 

TXRMOCOUPLE  2 

*• 

1.3012808 

>• 

*0.4281250 

THERMOCOUPLE  54 

M* 

1.0017774 

8* 

-0.45501 75 

T>«RMOCOUPL£  5 

H* 

1. 0022457 

8* 

*0.4744428 

T>OMOCOUPL£  37 

M* 

1.0015427 

8* 

-0.4770555 

T>«RMOCOUPLE  4 

Mt* 

1.0012808 

8- 

-0.4281250 

thermocouple  58 

M* 

1.0015207 

8* 

-0.5514854 

THERMOCOUPLE  5 

H- 

1 . 3015457 

B» 

-0.2541515 

TxeRKocautLe  st 

M* 

1.0024245 

1- 

thermocouple  4 

Mf* 

1.0015844 

8* 

-0.5557558 

THCRHOCOUPUE  <0 

M* 

1.0024214 

-0.7574074 

THERMOCOUPLE  7 

*m 

1.0017455 

8* 

-0.2480554 

TWRMOCOUPLE 

M* 

1.0017155 

B- 

THERMOCOUPLE  8 

1.0017471 

8* 

-0.5487588 

T>«HOCOU»L£  42 

M* 

1.0015044 

8* 

THERMOCOUPLE  5 

HI* 

1 .0024521 

8* 

-0.5058084 

THERMOCOUPLE  «j 

M* 

1.0015274 

8* 

-0.7550077 

THERMOCOUPLE  10 

HI* 

1.0015540 

8* 

-0.5740255 

T>«RHCCOUPLE  4* 

M* 

1.0024417 

8- 

-0.5244112 

TTERMOCOUPLE  11 

H- 

1.0017424 

8* 

-0.5452187 

THERMOCOUPLE  45 

M* 

0.5557850 

8* 

T5ERM0C0LTLE  IZ 

M- 

1.0024770 

8* 

-1.0715818 

THERMOCOUPLE  45 

M* 

1.0015827 

1* 

THERMOCOUPLE  15 

M* 

1.0024207 

8* 

-0.5058085 

THERMOCOUPLE  47 

H* 

1.0017757 

8* 

THERMOCOUPLE  14 

H* 

1.0015180 

8* 

-0.8755855 

THERMOCOUPLE  4# 

M* 

1.0050405 

B- 

THERMOCOUPLE  15 

H* 

1.0024844 

8* 

-0.8721151 

THERMOCOUPLE  45 

M* 

1.0025584 

8* 

-0.7045517 

THERMOCOUPLE  14 

Mt« 

1.0024144 

8* 

-0.8571475 

thermocouple  so 

M* 

1.0024405 

8* 

-0. 41 54777 

T>ERMOCOUPLE  17 

M* 

1.0021400 

B* 

-0.4405557 

THERMOCOUPLE  51 

M* 

1.0010144 

8* 

thermocouple  18 

H* 

1.0025555 

8* 

-0.8805175 

THERMOCOUPLE  52 

M* 

1.0015054 

8* 

thermocouple  If 

*• 

1.0018555 

8- 

-0.4084425 

THERMOCOUPLE  55 

M* 

1.0015405 

8- 

-0.4540537 

THERMOCOUPLE  20 

H* 

0.5514704 

8* 

2.2540512 

THERMOCOUPLE  54 

M* 

1.0013054 

8* 

rwe*Mocau»u£  21 

M* 

1.0018024 

8* 

-0.4140444 

THERMOCOUPLE  55 

M* 

1.0015588 

8* 

THERMOCQU PLE  22 

M» 

1.0015808 

8* 

-0.4554428 

THERMOCOUPLE  54 

M* 

1.0011254 

8* 

-0.547*222 

THERMOCOUPLE  25 

N* 

1.0011015 

8* 

-0.2445800 

T>«RMOCOUPLE  57 

M* 

1.0017183 

8* 

-0.5253417 

THERMOCOUPLE  24 

H* 

1.0015580 

8* 

-0.5554054 

thermocouple  St 

H* 

l .0-32502 

B» 

-0.8305081 

THERMOCOUPLE  25 

m* 

0.5584758 

8* 

0.5508555 

THERMOCOUPLE  5» 

M* 

l • 0024283 

-0. 7043454 

THERMOCOUPLE  24 

m** 

1.0022001 

8* 

-0.4275254 

thermocouple  ts 

M* 

1.0027475 

8* 

THERMOCOUPLE  2  7 

m* 

1.0018101 

8* 

-0.0544455 

THERMOCOUPLE  41 

M* 

1-0014108 

8* 

-0.4512557 

THERMOCOUPLE  28 

m* 

1.3026242 

8* 

-0.2547775 

THERMOCOUPLE  42 

M* 

1.0025404 

8* 

-0.7 730448 

THERMOCOUPLE  2f 

M- 

1.0010777 

8* 

-0.1212850 

THERMOCOUPLE  43 

M* 

1.0015532 

8* 

THERMOCOUPLE  50 

M* 

1.0045425 

8* 

-0.7454101 

THERMOCOUPLE  44 

M« 

1.0020733 

-0.4718745 

TVER  MOCQ  UP  L£  51 

M* 

1.0024007 

8* 

-0.4801757 

TXRMOCOUPLE  45 

M* 

1.0018148 

8- 

-0.5580448 

THERMOCOUPLE  52 

M* 

1.0010715 

8* 

-0.4208584 

THERMOCOUPLE  44 

M* 

1.0018148 

8* 

-0.5580648 

THERMOCOUPLE  55 

HI* 

0.5555580 

8* 

-0.1528710 

THERMOCOUPLE  47 

M* 

1.0018535 

8* 

-0.7208007 

THERMOCOUPLE  54 

M* 

1.0010715 

8* 

-0.4208584 

THERMOCOUPLE  4t 

M* 

1.0005155 

8* 

-0.5204545 

i 

\ 
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APPENDIX  B 


Heated 

TABULATED  DATA 

Cylinder  in  the  First  Row 

H  =  1 .0 

1 . 5m/s 

3  .  Om/s 

Downstream  Row  No. 

Theta 

Theta 

1 

1.0 

1  .0 

2 

.468 

.  524 

3 

.336 

.380 

4 

.251 

.311 

5 

.181 

.211 

Heated  cylinder  in 

the  first  row. 

H  =  4.6 

1 .  5m/s 

3  .  Om/  s 

Downstream  Row  No. 

Theta 

Theta 

I 

1.0 

1.0 

2 

.497 

.484 

3 

.365 

.380 

4 

.276 

.277 

5 

.  197 

.  185 

Heated  cylinder  in 

the  second  row . 

H  =  2.3 

1  .  5m  /  s 

3  .  Om/  s 

Downstream  Row  No. 

Theta 

Theta 

1 

1  .0 

1.0 

2 

.  501 

.514 

3 

.324 

.354 

4 

.250 

.269 

53 


TABULATED  DATA 


Row 

Ts.o(°C) 

T  (°C) 

0  ( w) 

h ( w/m^ 0  c  ) 

ReH 

Nu 

U  (m/s  ) 

No  . 

H  =  1 .0 

1 

98.5 

20.3 

3.03 

56.61 

394.2 

14.09 

1.5 

2 

98.9 

20.3 

3.05 

56.71 

394.2 

14.11 

1  .5 

3 

99.6 

20.5 

3.28 

61.16 

393.9 

15.21 

1.5 

4 

97.6 

20.0 

3.18 

60.34 

394.5 

15.02 

1.5 

6 

96.3 

20.1 

3.74 

73.87 

394.4 

18.39 

1  .  5 

1 

101.2 

20.8 

4.59 

87.27 

787.2 

21.69 

3.0 

2 

89.4 

20.7 

4.17 

93.26 

787.4 

23.19 

3.0 

3 

78.9 

20.5 

3.73 

98.52 

787.9 

24.51 

3.0 

4 

96.8 

20.2 

4.74 

95.26 

787.6 

23.71 

3.0 

6 

86.2 

20.3 

4.57 

107.69 

788.3 

26.80 

3.0 

H  =  2  .3 

1 

111.0 

20.0 

4.12 

45.77 

907.4 

26.21 

1.5 

2 

105.7 

20.1 

4.04 

47.93 

907.1 

27.44 

1.5 

3 

102.8 

20.1 

4.13 

51.03 

907.1 

29.22 

1.5 

6 

93.7 

20.1 

3.78 

52.61 

907 . 1 

30.12 

1.5 

1 

99.3 

20.2 

5.50 

73.19 

1813.7 

41.90 

3.0 

2 

110.1 

20.1 

5.99 

69.86 

1814.2 

39.99 

3.0 

3 

118.2 

20.2 

6.81 

73.20 

1813.7 

41.90 

3.0 

6 

71  .8 

20.1 

3.91 

80.03 

1814.2 

45.82 

3.0 

H  =  4 . 6 

1 

96.7 

20.0 

3.27 

42.74 

1814.7 

48.96 

1.5 

2 

97.3 

20.2 

3.27 

42.49 

1813.7 

48.65 

1.5 

3 

104.8 

20.3 

3.71 

44.20 

1813.2 

50.60 

1.5 

4 

100.7 

20.6 

3.00 

36.9 

1811.6 

42.22 

1.5 

6 

91  .5 

20.6 

2.92 

41.14 

1811.6 

47.06 

1.5 

1 

112.3 

19 .5 

5 . 60 

62.83 

3634.7 

72.04 

3.0 

2 

99.6 

19.8 

5.02 

65.70 

3631.6 

75.28 

3.0 

3 

95.6 

19.9 

4.80 

66.23 

3630.6 

75.87 

3.0 

4 

91  .2 

20.1 

3.92 

56.91 

3628.5 

65.17 

3.0 

6 

85.8 

20.1 

3.92 

61.93 

3628.5 

70.91 

3.0 
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